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Abstract 

We investigate ground state properties of a quasi-one-dimensional electron-lattice coupled model for quarter-filled 
molecular conductors. The effective one-dimensional extended Hubbard model coupled to adiabatic lattice degree of 
freedom is derived by the inter-chain mean-field approximation and solved by Lanczos exact diagonalization method. 
We find that the critical behavior among lattice tetramerized states with different charge-lattice ordered patterns is 
sensitively affected by the inter-chain Coulomb interaction, lattice anharmonicity, and intrinsic dimerization. This 
indicates a subtle balance between these states originating from strong correlation and reduced dimensionality. 
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1. Introduction 

Quasi-one-dimensional (QID) molecular conduc- 
tors exhibit a variety of phase transitions, where 
different kinds of symmetry-broken orderings in 
charge, spin, and lattice degrees of freedom are com- 
peting. Such competitions occur even when their 
non-interacting band structures are similar, sug- 
gesting a crucial role of interactions, i.e., electron- 
electron as well as electron-lattice interactions. 
Typical examples are the families of (DCNQI)2X 
and (TMTTF)2X, which have QID quarter-filled 
TT-bands in terms of electrons and holes, respec- 
tively [T]. They show different phase transitions 
depending on each material as well as applied pres- 
sure, such as charge ordering (CO), lattice dimer- 
ization, which drives the system effectively half- 
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filled resulting in a dimer Mott (DM) insulating 
state, and lattice tetramerization which is ascribed 
to a spin-Peierls (SP) transition |2I3I4I5I6I7| . This 
diversity originates from an interplay between the 
strong Coulomb interaction among electrons and 
the reduced dimensionality giving rise to lattice 
instabilities. 

Theoretical aspects of such electron-lattice cou- 
pled phase transitions have been extensively stud- 
ied based on ID and QID extended Hubbard 
models with electron-lattice couplings at quarter- 
filling [5]. Different states have been reproduced, 
at the ground state [911011111^ as well as at finite 
temperatures |13I14I15I16I17] . In refs. [TTTT^ . the 
ground state phase diagrams of ID models were ob- 
tained, where two different kinds of SP states were 
found to be competing with each other. One is the 
coexisting state with CO (CO-I-SP) and the other is 
that with lattice dimerization, i.e., DM (DM-fSP). 
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The phase boundary between these two states is 
of first order. On the other hand, the authors have 
recently investigated finite-temperature properties 
of a more general model including the inter-chain 
Coulomb interaction, a lattice anharmonicity, and 
the intrinsic dimerization, and found that they play 
crucial roles in the critical behavior among different 
phases, and are relevant to the experimental sys- 
tems (17j . In this paper, motivated by the results, 
we investigate the effects of such additional con- 
tributions on the ground-state phase competitions 
and coexistences, and how the criticality among 
different SP states mentioned above are affected. 

2. Formulation 

We investigate the QID quarter-filled model pre- 
viously studied in refs. |15I17| at finite temperatures, 
whose Hamiltonian is given by H — J2j ^id + '^-^ ' 
where Ti^p and Ti.± are the intra-chain part of the j- 
th chain and the inter-chain part, respectively. The 
former is given by the extended Peierls-Hubbard 
Hamiltonian, 

X (4j\<tC»+1 J,^ + h.c.) 

i i 

(1) 

i i 

where the notations are referred to ref. [17j . Here we 
consider the uniaxial lattice distortions Uij along 
the ID chain and neglect their quantum fluctuations 
as in the previous studies |llll2j . The term with 
-fCpj represents an anharmonicity in the lattice dis- 
tortion. The inter-chain part is given by 

= Vj_ E E '^i.l'^i.k^ (2) 

where V± denotes the inter-chain Coulomb interac- 
tion and the summation for (j, fc) runs over nearest- 
neighbor chains. In the following calculations, we 
choose the on-site and nearest-neighbor Coulomb in- 
teractions as U — 6 and V — 2.5, respectively, and 
the elastic constant as Kp = 0.75, in the energy unit 
oft [TT]. 



(a) CO+SP 




Fig. 1. Schematic view of the interchain configuration 
between neighboring chains we consider in this work, for 
(a) coexistence of charge order and spin-Peierls (CO+SP) 
states and (b) coexistence of dimerized Mott insulator and 
spin-Peierls (DM-I-SP) states. Size of the sites and thickness 
of the bonds represent the charge density and the lattice 
distortion, respectively. The figure is drawn for the two-di- 
mensional case (z = 2) where chain j is coupled to chains 
j — 1 and j + 1. 

We derive an effective ID model by the inter-chain 
mean- field approximation as in ref. [17j . The mean- 
field form of H± is taken as 

i 

- (3) 

where z denotes the number of nearest-neighbor 
chains and the chain index j is dropped hereafter. 
Eq. ([3]) is derived by assuming the relative phase 
between neighboring chains as drawn in Fig. 1. We 
treat Ui adiabatically, then static lattice distortions 
are determined so as to minimize the ground-state 
energy under the constraint — 0- 

We obtain the ground state of the effective ID 
model TCiB + TC^^ by exact diagonalization using 
the Lanczos technique, and determine (rii) and Ui 
self-consistently. We investigate symmetry break- 
ing with a four-sites period along the chain, taking 
into account all possibility of the CO, DM, and 
their variants with further tetramerizations. The 
self-consistent equations then involves the charge 
densities {ni) and the lattice distortions ui for 
I — 1,2,3,4. In the following, we show results for 
the system size of L = 12 under the anti-periodic 
boundary condition [18 . 

The obtained self-consistent solutions are param- 
eterized as (see Fig. 1), 

(ni) = n - nco ~ ^it(co+SP) ~ ^it(DM+SP)i 
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(n2)=n + rico +"-t(DM+SP), 
(ng) = fi - nco + "•t(co+SP) + "t(DM+SP), 
(n4)=n + nco -«t(DM+SP)7 

til = +Md - '"t(CO+SP), 

U2 = -Md + '"t(C0+SP) + "t(DM+SP)j 

"3 = +""£1 + Ut(CO+SP), 

U4 = -Md - Ut(CO+SP) - '"t(DM+SP)7 (4) 

where n = 1/2 is the average density, nco repre- 
sents the alternating charge disproportionation for 
the CO state, and Ud the lattice dimerization for 
the DM state, respectively. Lattice tetramerizations 
described by Ut(co+SP) and Mt(DM+SP) characterize 
the CO+SP and the DM+SP states, respectively, 
while four-fold charge disproportionation concomi- 
tantly appears as nt(co-i-SP) a^nd ntpM+sp) in each 
of these states. 

3. Results 

First, we show the results in the absence of intrin- 
sic dimerization, i.e., Sd = 0. Fig. 2 shows the vari- 
ation of the order parameters as a function of zVj_ . 
For zV± < 0.88 the ground state is the DM-fSP 
state, where {ud, Mt(DM+SP), «t(DM+SP)} are finite 
while nco = «t(co+SP) = '«t(co+SP) = 0. On the 
other hand, the CO-I-SP state appears for zV± > 
0.93 where, vise versa, {nco, '^t(co+SP) , ■"t(co+SP)} 
are finite and Ud = Ut(DM+SP) = «t(DM+SP) = 0. In 
the narrow intermediate region, 0.88 < zV± < 0.93, 
all of the six parameters become finite. This indi- 
cates that a further coexistence of the DM+SP and 
CO-I-SP states is stabilized, which was not found 
in previous studies on the ground state properties 
for similar models [ll|12j . The two phase transition 
points, along DM-fSP the multiple coexistent 
state CO+SP, are of second order. 

Next, for the case of anharmonic lattice distor- 
tions with Kp^ = 0.75, the results are shown in 
Fig. 3. Here, we find different critical behavior be- 
tween the DM+SP and CO+SP states; they are 
bordered by a direct first order phase transition at 
zV± = 0.74. This is interesting in a sense that, at fi- 
nite temperatures above the SP transition tempera- 
tures, our recent results show that the anharmonic- 
ity contributes to rather enlarge the coexistent re- 
gion of the paramagnetic CO and DM states [T7]. 

Finally, the results in the presence of intrinsic 
dimerization Sd are shown in Fig. 4; here, Ud is al- 
ways finite due to the symmetry breaking from the 




Fig. 2. Order parameters as a function of zV± for (5^ = 
for the harmonic case with Kp^ = 0. 

outset in the model. For zV± < 1.0 the ground 
state is in the SP(+iDM) phase with {ut(i3M_|_gp), 
"t(DM+SP)} being finite and nco = "^t(co+SP) = 
''^t(CO+SP) = 0. Here iDM refers to the intrinsic 
dimerization to distinguish with the spontaneous 
dimerization in the model with Sd — 0, although the 
actual pattern is the same as in the DM+SP state. 
As for zV± > 1.0, all the six order parameters be- 
come finite; there occurs a second order phase tran- 
sition into the CO+SP (+iDM) state 17]. As a re- 
sult of the intrinsic dimerization leading to Ud ^ 0, 
when the CO+SP component arises the order pa- 
rameters for the SP(+iDM) state always remains, 
then the charge-lattice pattern is the same as in the 
multiple coexistence in Fig [51 
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Fig. 3. Order parameters as a function of zV± for 
for the anharmonic case with Kp = 0.75. 
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Fig. 4. Order parameters as a function of zV± with intrinsic 
dimerization = 0.02 for the harmonic case with Kp^ = 0. 

4. Discussions 

Here we compare our results with those in the 
previous studies. In ref. ^11 , the model with = 
and zV± = was studied. There, a first order 
phase transition between the DM+SP and CO+SP 
states was found, while the multiple coexistence was 
not. This may be ascribed to the difference between 
the models, i.e., whether to include zV± or not. 
On the other hand, a model including Holstein-type 
electron-lattice coupling was considered in ref. [12] , 
which also indicated the first order phase transition. 
The Holstein coupling under the adiabatic approx- 
imation leads to a term with a similar but slightly 
different form compared with the inter-chain mean- 
field contribution considered in the present work. 
This may be the cause of the difference in the re- 
sults. All these comparisons suggest that the inter- 
chain Coulomb interaction as well as lattice anhar- 
monicity affects the critical behavior in the electron- 
lattice coupled phase transitions. A comprehensive 
study including the system size extrapolation is left 
for future study. 

5. Summary 

In summary, we have investigated the phase 
competition and coexistence between different spin- 
Peierls phases observed in quasi-one-dimensional 
molecular conductors. We found different sequences 
of phase transitions at the ground state, which re- 
flect the slight difference in the model. This suggests 
a subtle balance between charge-lattice coupled 



symmetry breakings arising from strong correlation 
and reduced dimensionality. 
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